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Abstract. We study discrete Schrodinger operators with compactly supported 
potentials on Z"*. Constructing spectral representations and representing S- 
matrices by the generalized eigenfunctions, we show that the potential is 
uniquely reconstructed from the S-matrix of all energies. We also study the 
spectral shift function 5(A) for the trace class potentials, and estimate the 
discrete spectrum in terms of the moments of 5(A) and the potential. 

1. Introduction 

In this paper, we consider trace formulas and inverse scattering problems for 
Schrodinger operators on the square lattice Z"^ with d> 2. We restrict ourselves to 
the case of compactly supported, or trace class, potentials. Our aim is two-fold : 
the recontsruction of the potential from the scattering matrix, and the computation 
of trace formula using the spectral shift function. We begin with the forward 
problem. We shall construct the generalized Fourier transform and represent the 
S-matrix by generalized eigenfunctions. We then show that given the S-matrix 
S{X) for all energies, one can uniquely reconstruct the potential (Theorem 4.4). We 
next compute the asymptotic expansion of the perturbation determinant associated 
with the discrete Hamiltonian. By virtue of Krein's spectral shift function, one 
can compute the moments of logdetiS(A). As a by-product, one can estimate the 
discrete spectrum using these moments (Theorem 6.4). 

In the continuous model, the first mathematical result on the inverse scattering 
for multi-dimensional Schrodinger operators was that of Faddeev |13| : the recon- 
struction of the potential from the high-energy behavior of the scattering matrix 
using the Born approximation. In this paper, we give its discrete analogue. Instead 
of high-energy, we consider the analytic continuation of the S-matrix with respect 
to the energy parameter and use the complex Born approximation. The analytic 
property of the resolvent of the discrete Laplacian is more complicated than the 
continuous case, which requires harder analysis in studying the inverse scattering 
problem. In the continuous case, Faddeev proposed a multi-dimensional analogue of 
the Gel'fand-Levitan theory using new Green's function of the Helmholtz equation 
([15], [16]). Faddeev's Green function was rediscovered and developed in 1980's 
by Sylvester- Uhlmann [33], Nachman [27], Khenkin-Novikov [25] (see the survey 
article [20] )■ Reconstruction of the potential from the S-matrix of a fixed energy 
is one of the novelties of this idea. In the forthcoming paper, we shall study the 
inverse scattering from one fixed energy in the discrete model. 

Our next purpose is the trace formula. It is well-known that the scattering 
matrix, Krein's spectral shift function and the perturbation determinant for a pair 
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of self-adjoint operators Hq and H ~ Hq + V are mutually related. We shall write 
down the first 5 moments of the spectral shift function for our discrete model in 
terms of the traces of the potential V. If the potential admits a definite sign, we 
can obtain estimates of the discrete spectrum by the traces of V. 

The computation of trace constitutes the central part of the study of spectral 
theory, since it provides quantitative information of the operator in question, hence 
serves as a clue to the inverse problem. In the continuous case, the trace formula 
was first obtained by Buslaev-Faddeev [TT] in the one-dimensional problem and 
by Buslaev [9], |10| in the three-dimensional problem. Since then, an abundance 
of articles have been devoted to this subject, see e.g. [H] and [H], [5D] . 
|17) . Gesztesy-Holden-Simon-Zhao [15] computed the trace Tr (L — L^), where 
L = —A + V is a Schrodinger operator for the continuous case, and is L with 
Dirichlet condition on a subset A C R'^. Shirai [32] studied this problem on a 
graph. Karachalios [23] and Rosenblum-Solomjak [31] computed the Cwikel-Lieb- 
Rosenblum type bound for the discrete Schrodinger operator. The well-known 
Efhmov effect has a different property in the case of the discrete model. See e.g. 
Albeverio, Dell Antonio and Lakaev [2]. See also [3]. 

In §2, we shall prove the limiting absorption principle with the aid of Mourre's 
commutator theory |26j. We then derive the spectral representation in §3, and rep- 
resent the S-matrix by generalized eigenfunctions. The inverse scattering problem 
is solved in §4, and §6 is devoted to the trace formula. 

The essential spectrum of our discrete Schrodinger operator H — Hq + V on Z"^ 
fills the interval [0,d], and the set [0,d] H Z is that of critical points, since Hq is 
unitarily equivalent to the operator of multiplication by {d — X]j'=i cosxj)/2 on the 
torus R"*/ (27rZ)'^. In fact, the resolvent estimates in §2 are proved outside [0, d] D Z. 
The behavior of the free resolvent {Hq — z)~^ near the critical points depends on 
the dimension d. This is itself interesting and is studied in §5 (Lemmas 5.3, 5.4 and 
5.5), although we do not use it in this paper. 

The notation used in this paper is standard. For two Banach spaces X and Y, 
'B{X; Y) denotes the set of all bounded operators from X to Y. For a self-adjoint 
operator A, cr{A), (Jp{A), ad{A), ae{A), aadA) and p{A) denote its spectrum, 
point spectrum (= the set of all eigenvalues), discrete spectrum, essential spectrum, 
absolutely continuous spectrum and resolvent set, respectively. For a trace class 
operator K, Tr (K) denotes the trace of K. 

2. Schrodinger operators on the lattice 

2.1. Discrete Schrodinger operator. Let Z"^ = {n = (ni, - - ,71^); G Z}, 
and ei = (1, 0, • • • , 0), • • • , — (0, ••• ,0,1) be the standard bases of Z"^. The 
Schrodinger operator H on Z*^ is defined by 

H = Ha + V, 

where for /= {/(n)}^^^. G PiZ'') and n e Z'' 

{Haf} (n) ^Pin)-\J2 + + " ^^)} ' 

iVf)in) = Vin)fin). 
Until the end of §4, we impose the following assumption on V: 
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(A-1) V is real-valued, and V{n) — except for a finite number of n. 
Define the l-dimensional projection P{n) by 

(P(n)7) (m) = 5n^f{m). 

Then V is rewritten as 

V ^Y.^{n)P{n). 

n 

Let us introduce the shift operator 

{S,f){n) - f{n + e,), f){n) = f{n - e,), 

and the position operator 

{N,f){n)=nJ{n). 
A direct computation yields the following lemma. 

Lemma 2.1. Sj is unitary, Nj is self-adjoint with its natural domain and 
Letting N — {Ni, ■ ■ ■ , Nd), H is rewritten as 

i=i 

The spectral properties of H is easier to describe by passing to its unitary trans- 
formation by the Fourier series. Let 

T'' = RV(27rZ)'' = [-TT,Trf 

be the flat torus and U the unitary operator from P{Z'^) to L^(T'') defined by 

{Uf){x) = {2n)-''/' /(n)e— . 

neZ'' 

The shift operator and the position operator are rewritten as 

Sj-.^U U* = e"^ , Nj-.^U Nj U* = id/dxj . 

Letting 

Ho^UHqU*, V=UVU*, 

we have 

(2.1) H ■=UHU* = Ho + V, 

1 

(2.2) Ha = -(d-Ycosxj'j =: h{x), 

iVf){x) - (2^)-'^/2 / V{x - y)f{y)dy, 
V{x) = (2^)^^/2 T/(n)e-*"-". 

In fact, this is a special case of the Friedrichs model (see e.g. [2]). The following 
theorem is easily proven by (|2.ip and Weyl's theorem. 
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Theorem 2.2. (1) a{Ho) = aadHo) = [0,d]. 
(2) ae{H) = [{),dl C R\ [0,d]. 

2.2. Sobolev and Besov spaces. Wc put N = {Ni, • • • , Nd), and let iV^ be the 

self-adjont operator defined by 

d 

= Y^ = -A, on T'^, 

where A denotes the Laplacian with periodic boundary condition. 

We put 

|iV| = \/7V2 = 

We introduce the norm 

\\u\\s^\\{l + Ny/''u\l seR, 

II • II being the norm on L'^{T'^), and let W be the completion of D{\N\''), the 
domain of \NY , with respect to the norm ||m||s : 

W = {ueV'{T'') - ||m||.<oo}, 

where V'iT^) denotes the space of distribution on T''. Put U = = L'^ {T^) . 

For a self-adjoint operator T, let x{a < T < b) denote the operator xi{T), 
where X/(A) is the characteristic function of the interval / = [a, 6). The operators 
x{T < a) and x{T > b) are defined similarly. Using the sequence {rj}J^Q with 
r_i =0, rj ~ 2^ (j > 0), we define the Besov space B by 

oo 

B={/eH;||/||B = ^r]/2||x(r,-i < \N\ < rj)f\\ < ^} . 

Its dual space B* is the completion of H by the following norm 

||w||g. = supr7^/2||x(rj-i < |A^| < rj)u\\. 

The following Lemmas 2.3 and 2.4 are proved in the same way as in the contin- 
uous case (|T]). We omit the proof. 

Lemma 2.3. There exists a constant C > such that 

/I \ 

C\\u\\b^ < sup -||x(|iV| < R)uf < C\\u\\b'. 

\R>1 ^ J 

Therefore, in the following, we use 

/I 

Mb- - (^sup-||x(|iV| <i?)u||2j 

as the norm on B* . 

Lemma 2.4. For s > 1/2, the following inclusion relations hold : 
cBc H^/^ cHc n-^'"^ cB* c . 

We also put 11 = /^(Z''), and define W , B, B* by replacing N by N. Note that 
— U*H'^ and so on. In particular, Parseval's formula implies that 

M'u^ = M%= ^(l + H2)^Rn)p, 
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uWa, = sup \u{n)f, 



|n|<-R 



u(n) being the Fourier coefficient of u(x). 

2.3. Mourre estimate. Let — Dsyo'H'', and define a symmetric operator A 
with domain 71°° by 

(2.3) A = ^J2 [{iS^r - S,)N, + ^1±^^ . 

Then A is essentially self-adjoint. In fact, letting M = 1 + N'^, we can find a 
constant C > such that 

(2.4) \\Au\\<C\\Mu\\, \{Au,Mu)-{Mu,Au)\<C\\M^^'^uf, yueH°°. 

Nelson's commutator theorem ([2^, p. 193) then implies the result. 
By Lemma 2.1, (|2.3I) is rewritten as 



Let us note that in [8^, — (^-^j^j is used as A. Our choice of A 

comes from the following reasoning. Let h{x) be defined by (|2.2p . Passing to the 
Fourier series, we have 



d „ 
a 



A = UAU* = (2sinxj- h cosx^j = 2i{Vxh -Vx + ^x ■ (V^^)) 



This is an analogue of the generator of dilation group on R"*. We then have 

d 

i[Ho,A]=4\Vxh\^ ^Y.'^sinx,)^. 

Let Eho{)^) and Eh{X) be the spectral decompositions of Ho and H, respectively. 

Lemma 2.5. Let A G (0, d) \ Z. Then there exist constants (5, C > and a compact 
operator K such that 

Eh{I)[H, iA]Eh{I) > CEh{I) +K, I^{X~d,X + S). 

Proof. For A e (0, d) \ Z, let 

Mx = {x e h{x) = A}. 

If V/i(x) = 0, then cosa;j — ±1, and h{x) e Z. Therefore, the assumption A ^ Z 
implies that Vh{x) 7^ on Mx, hence M;)j is a real analytic manifold. We put 

Co(A) = inf |V/i(x)p. 

Then for any small e > 0, there exists 5 > such that 

\Wh{x)f >Co{X)-e on /i-i([A - (5, A + (5]). 
We have, therefore, 

EHoiI)[Ho,iA]EHoiI) > (CoiX) ' e)EH„iI), / = (A - <5, A + (5). 
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Since ^ is a compact operator, so are V and Eh {I) — Ehq{I)- This proves the 
lemma. □ 

Let R{z) = {H — z)^^ be the resolvent of H . 

Theorem 2.6. (1) (jp{H) n ((0, d) \ Z) is discrete and finite multiplicities with 
possible accumulation points in Z. 

(2) Let s > 1/2 and A £ (0,d) \ (Z U (Tp(i/)). Then, thers exists a norm limit 
R{X±iO) := lim^^o R{X ± ie) e B{W-H-'). Moreover, we have 

(2.5) SUp||i?(A±iO)||B(6;6-) <f», 

AGJ 

for any compact interval J in (0,c?) \ (Z U (Jp{H)'^ . The mapping (0, rf) \ (Z U 
{H)) 3 B.{X±iQ) is norm continuous in 'B{'H'^;'H and weakly continuous 
mB{B;B*). 

(3) H has no singular continuous spectrum. 

This theorem follows from the well-known Mourre theory. We shall give here a 
brief explanation. First we introduce an abstract Besov space. We define 

oo 

BA^{feH; \\f\\B. = ^rj/'||x(r,-i < |^| < r,)/|| < oo}, 

where TL = L'^{T'^). Its dual space Ba* is the completion of TL by the norm 

\\u\\ba' ^suprJ^^^Wxirj^i < \A\ < rj)u\\. 
j 

The abstract theory of Mourre based on Lemma 2.5 then yields 

Lemma 2.7. Let J be as in Theorem 2.6 (2). Then there exists a constant C > 
such that 

sup \\{H-zr'f\\B,-<C\\f\\B,, yfeBA. 

For the proof of the lemma, see j4], j8], [22]. Therefore, to prove Theorem 2.6, 
we have only to replace Ba by B using the following lemma. 

Lemma 2.8. There is a constant C > such that 

ll/b. <C||/||b, V/eS. 

Proof. For t G R, let (t) = (1 + 1^)^/^. By definitions of A and N we have 

{A){N)-^ e B{n;n). 

For feB,we put /, = x{rj^i < \N\ < r,)f. Then 

||(A>/,|| <q|(7v)/,|| <ar,||/,|i, 

which implies 

||x(r,-i < 1^1 < rfe)/,|| = ||x(rfe-i < |^| < rk){A)-^ {A) f,\\ 
<Cr-'\\{A)f,\\<Cr-'r,\\f,\\. 

Then we have 

Y.^l^'\\xirk-i < \A\<r,)f,\\<Cj2r-'/'r,\\f,\\<Crf\\f,\\, 

k>j k>o 
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Y.^l''Mrk-i < \A\ < r,)f,\\ <J2ri^'m\ < C'rfWm- 
k<j k<j 

We have, therefore, 

<Cr]/^ll/,ll, J=0,l,2,.... 
Summing up these inequahties with respect to j, we obtain the lemma. □ 

3. Spectral representations and S-matrices 

3.1. Spectral representation on the torus. For t E {0,d) \ Z, let dMt be the 

mesure on Mt induced from dx. By taking t = h{x) as a new variable, one can 
show that for / e C(T'*) supported in {a; G T'* ; h{x) ^ Z} 

(3.1) 



For /,g e L^j-rjid-)^ ^^^^ 



d h{x) — z 

Therefore, if /, g € C^{T'^) and A G (0, d) \ Z, 

(3.2) (i?o(A±zO)/,5) = ±*^ / fg-^+p.y.f ,.{^ . dx. 

Let L'^{Mx) be the Hilbert space equipped with the inner product 
(3-3) (v3,V')l2(Ma) = / (pi^^^^. 

We define 

(3.4) J^o(A)/ = / , 

where the right-hand means the trace on , i.e. the restriction to, M\. Then we 
have by ([^ 

Lemma 3.1. For \ e (0, d) \ Z, and /, g E (T'^), 

±-,{{R,{X + ^0)-R,{X-^Oj)f,g) = {M^}f,M>^)9)L2^Mxy 
By (12. 5p . this lemma implies 

(3.5) MX) EB{B-L\Ah)). 

Moreover, 

(3.6) {f,g)L2^T'')^ [ {MX)f.MX)9)LHMx)d\ f,geB. 

Jo 

The adjoint operator Tq{X)* is defined by 

(•^o(A)/, 0)l2(m^) = (/, J"o(A)»L2(xd). 

By (1331), Jo(A)* e B(L2(M^):S*), and by To{X){Ho - A) = 0. Hence we 

have 

(i/o - A)J-o(A)* =0. 
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In view of we can identify L^(T'^) with the space of L^-functions /(A) over 
(0, d) with respect to the measure dX such that for a.e. A £ (0, d), /(A) takes values 
in L'^iMx). We denote this space by £^((0, d); L^{Mx); dX). 

We put (Jo/) (A) = J^o(A)/ for f e B. The following Theorem 3.2 is essentially 
a reinterpretation of the identification L'^{T^) ~ L^{{0,d); L^{AI\);dX). However, 
we give a functional analytic proof for the later convenience. 

Theorem 3.2. (1) is uniquely extended to a unitary operator 

To ■■ L'iT'-) ^ l2((0, d);L\Mx);dX). 

(2) diagonalizes Hq : 

(J-oi/o/)(A) = A(Jo/)(A), yfeL^T''). 

(3) For any compact interval I C (0, d) \ Z, 

J-o(A)*5(A)rfA e L^T"), Vg e L\{0, d);L^{Ahy, dX). 
Moreover, for In = U^=i(i - 1 + 1/^^, j - ^/N), the inversion formula holds: 
/= lim / J-o(A)*(J-o/)(A)dA, yfeL^T"), 

J -In 

where the limit is taken in the norm of L^{T^). 

Proof. By (|3.6p . J-q is uniquely extended to an isometric operator from L'^{T'^) 
to L^((0, d); L^(AfA); dA). To show that it is onto, we have only to note that the 
range of is dense. For f £ B, we have J^o{X){Ho — A)/ = by definition, which 
proves (2). To show (3), we first note that for a compact interval / C (0,d) \ Z, 
Jj To{X)* g{X)dX e B*. We use ( , ) to denote the inner product of L^(T'*) as weU 
as the coupling of B and B* . Then we have 



J-o(A)*5(A)rfA,/j = J^iM^rg{X),f)dX 

^ J^ig{X),ToWf)mM.^dX, feB. 
Therefore 

(^^ J-o(A)*5(A)dA,/) < . llJ-o/ll = 11.911 • |1/|1. 
By Riesz' theorem, we then have 

J^ToW*g{X)dX e L^{T^), \\ J^ToW*giX)dX\\ < \\g\\. 
Therefore for any compact interval J C {0,d) \ Z, 

J-o(A)*(J-o/)(A)dA|| < ||J-oX,7(i/o)/ll = \\xj{Ho)f\\, 



where xj is the characteristic function of J. One can then show the existence of 
the strong limit 

lim / J-o(A)*(J-o/)(A)dA-:/. 

J 1m 
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By Parseval's formula, we have for any h S L^(T'*) 



{f,h) = lim / (J-o/(A),J-o(A)/i)dA 

J In 

= i^of,:Foh) = if,h), 
which imphes f = f. □ 

Next let us construct the spectral representation for H. We put 

^■(^^(A) = J'o(A)(l- V'i?(A±iO)), Ae {Q,d)\{ZU(jp{H)). 
Then by (|231l 

Lemma 3.3. For A e (0, d) \ (Z U crp{H)), and f,g eB 

((i?(A + zO)-i?(A-zO))/,.g) - (j-W(A)/,-F(±)(A).9) 

Proof. We put 

Hi^Ho, H2 = H, R,{z)^{H,-z)-\ 



1 

Then we have, by the resolvent equation 



E'^{X}^—{R,{\ + iO)-R,{\-iO)). 



^ {Rk{X + te) - RkiX - lej) 



:Gjk{X ±^e)*— (i?j(A + ie) - R,{X - ie)) G,k{X ± ie). 



1 

Letting e — )■ 0, we have for f,gEB 

(3.7) {El{X)f,g) = (£;;(A)G,fc(A±*0)/,G,fe(A±zO)5). 

Let j = 1, k = 2. Since J^(*'(A) ~ J'o(A)Gi2(A ± iO), the lemma follows if we 
replace f,g in Lemma 3.1 by Gjk{X±iO)f, Gjk{X ± iO)g. □ 

We define the operator J-(±) by (J"(±)/)(A) = for f E B. Let ■Hac(H) 

be the absolutely continuous subspace for H. 

Theorem 3.4. (1) J^*^*^ is uniquely extended to a partial isometry with initial set 
TLaciH) and final set L'^{{0,d);L^{Mx);dX). Moreover it diagonalizes H : 

(3.8) (j-Wi//)(A)=A(j-(±)/)(A), V/eX2(T'^). 

(2) The following inversion formula holds: 

(3.9) / = s-lim / J-(±'(A)*(-^'*V)(A)dA, VfeHac{H), 

where In is a union of compact intervals C (0,0?) \ (ZU o'p(if)) such that In — >■ 
(0, d) \ (Z U a-p{H)) as N ->-oo. 

(3) J-^^\X)* G 'R{L'^{M\)]B*) is an eigenoperator for H in the sense that 

[H -X)F^^\xy(t)^Q, y(j)E L^{Mx). 
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(4) The wave operators 

(3.10) W^^'^ = s - lim e''"e~''"" 

exist and are complete, i.e. the range of W'-^^ is equal to HaciH)- Moreover, 

(3.11) W^(±) = (-F(±))*J-o. 

Proof. The proof of (1), (2), (3) is the same as that for Theorem 3.2 except for 
the surjectivity of J^^'^\ Since V is trace class, the existence and completeness of 
wave operators p.lOp can be proven by Rosenblum-Kato theory (see [23], p. 542). 
The relation p. lip is also well-known, and we omit the proof (see e.g. [21] )• We 
then have J^'^^ = J^o (VF*^^') . The completeness of VF^*-* implies the surjectivity 

Ofj-^i). □ 

3.2. Spectral representation on the lattice. Theorem 3.4 is transferred on the 
lattice space by U. Wc put J-o{)^), J-'^'^^X), J-q and by 

J^o(A) - J-o(A)W, .f(±)(A) - J-(±)(A)Z^, 

We also define 

Theorem 3.5. (1) J-^^^ is uniquely extended to a partial isometry with intial set 
'H-ac{H) and final set L'^{T^). Moreover it diagonalizes H : 

(3.12) (.F(±)i//)(A) = A(^(±)/)(A). 

(2) The following inversion formula holds: 

(3.13) / = s-lim/ ^(±)(A)* (A)dA, yfeHadH), 

where In is a union of compact intervals C (0,d) \ U (7p{H)j such that In — >■ 
(0, d)\(Z(J ap{H)j as N -^oo. 

(3) J-^^\X) € B(L^(A/a) ; 13*) is an eigenoperator for H in the sense that 

(4) The wave operators 

=s-lime'*^e-"^° 

t— >±oo 

exist and are complete. Moreover, 
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3.3. Generalized eigenvector. It is customary to define the distribution d{h{x) — 
X) e 2?'(T'^) by 

f{x)S{h{x) - X)dx = / fix) J^' f e C{T''). 
T'i Jmx Nxh[x)\ 

We then see that J^o{\)* defines a distribution on T'^ by the following formula 

J'o{\)*(l) = (l){x)S{h{x) ~ X). 

The right-hand side makes sense when, for example, (f> e C°°{M\) and is extended 
to a C°°-function near M\, which is denoted by (f>{x) again. The Fourier coefficients 
of J-o{X)*(j) are then computed as 

(3.14) (27r)-'^/2 f e^--^S{h{x) - X)cl,{x)dx = {2n)-^^' [ 

We look for a parametrization of Mx suitable for the computation in the next 
section. Let us note 

-(d-^cosx,) =^sin2 (^), 
i=i i=i 

which suggests that the variables y — {yi, - ■ ■ ,yd) & [—1,1]'*: 

yj = sin , Xj — 2 arcsin yj 

are convenient to describe Hq. In fact, the map x — >■ ?/ is a diffeomorphism between 
two tori: 

RV(27rZ)'' = [-TT, nf^x-^ye [-1, 1]'* = RV(2Z)'*. 

Conseqently, 

^d-l 



(3.15) x{Vxe) = (2 arcsin(\/A6li), • • • , 2 arcsin(VA6'<j)) , 

gives a parameter representation of Mx. Passing to the polar coordinates y = y/X9, 
we also have 

(3.16) dx = J{y)dy = '^-^^^itl^^dXde , 
which implies 

(3.17) ^ ^ iVXt^JiVXO)^^^ 
^ ' |V,/i(a;)| 2 

We define i/^f") (n. A, 6*) by 

V^(o)(n, A, 9) (2vr)-'*/2i^^^^ e"-^(^«) .1(^X9) 

- xiVX9) 



n^licos {xj{VX9)/2)' 

where xiv) is the characteristic function of [—1, l]"*. By p.l4p and (I3.17p . we have 
for e L^{Mx) 

(3.19) {^oW*(f>){n)= [ ^^°\n,X,9)(f>{x{VX9))d9. 
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One can also see that if / is compactly supported 

(3.20) {ToWf)ixiVX9)) - (27r)-'^/2 ^ e-™--(^^)/(n). 

3.4. Scattering matrix. The scattering operator S is defined by 

S = {W+)*W-. 
We conjugate it by the spectral representation. Let 

S = J-Q S (J^o) , 

which is unitary on L^{{0, d); L^{M\); dX). Since S commutes with Hq, S is written 
as a direct integral 

S = [ ®SiX)dX. 

J{0,d) 

The S-matrix, 5(A), is unitary on L'^{M\) and has the following representation. 
Theorem 3.6. Let X G (0, d) \ (Z U ap{H)) . Then S{X) is written as 

5(A) = l-27ra(A), 

(3.21) A{X) ^MX)VMXy -M>^)VR{>^ + iO)VM>^y■ 
S'mce the proof is well-known, we omit it (see e.g. [21] )• 

By p.l8|) and (I3.20[) . the first term of the right-hand side of (|3.2ip has an integral 
kernel 

(3.22) (27r)-'^^J(VA0)J(VA0') ^ ^-^n■{x{V^e)-x(V^e'))y^^y 

The 2nd term of the right-hand side of p.2ip has the following kernel 
(3.23) 

„ (2^)-^/2 i^Y''^ j^^ff-^ e-"'-"(^^)y(n) (i?(A + iO)F^(°)(A,0')) H- 

4. Inverse scattering 

In this section we prove that the potential V is uniquely reconstructed from the 
scattering matrix for all energies. We first consider the analytic continuation of 
x{y/Xe) defined by (l3l^ . 

Lemma 4.1. Let t be a constant such that —1 < t < 1, r 7^ 0. Then f{z,T) = 
2arcsin(zr) is analytic with respect to z Cz C\ ^(— c», — 1/|t|] U [1/|t|, cx))^ . More- 
over, letting e{T) — t/\t\, we have as N ^ 00, 

Ref{N + t,T) = e{T)(n-^ + 0{^)), mod 27rZ 

Im f{N + z, r) = e(T) (2 log N + \og{4T^) + 0(^)) • 



DISCRETE SCHRODINGER OPERATORS 



13 



Proof. Wc take the branch of arcsm(z) = u{z) + iv{z) so that it is single- valued 
analytic on C \ ((-oo, -1] U [l,oo)), and < u{s) < 7r/2 when < s < 1. Then 
we have 

(4.1) sin(u) cosh(w) = Re z, cos(w) sinh(w) = Imz. 
Let us note that 

(4.2) cos(u) > and ± sinh(w) > 0, if ±Im2>0. 

In fact, by the 2nd equation of (|4.1|) . cos(m(2:)) and sinh(w(z)) do not vanish if 
Imz 7^ 0, and cos(m(s)) > 0, v{s) — when s G (0,1). So, cos(u(z)) > when 
Imz 0, and again by the 2nd equation of (j4.ip . sinh(u(z)) and Im z > have the 
same sign. 

Let UN = u{{N + i)T),VN — v{{N + i)T). Then by (|4.1|) we have sinh(t;Ar) — 
t/cos{un). Plugging this with 

(1 - COS^(uAr)) (1 + sinh^(t;Ar)) = TVV^ 

and letting tj^ = cos^(uAr), we get the equation 

t% + (iVV^ + - l)tN - = 0. 

Since > 0, by solving this equation, we have 

tN = N-^ + 0{N-^). 

Since cos{un) > 0, we have 

(4.3) cos{un) = + 0{N-^). 
This, combined with (|4.ip for z = {N + i)T, then yields 

(4.4) smh{vN) =tN + 0{N-^). 

If r > 0, then vn > 0, and by (|4T|) with z = {N + i)T, sin(ujv) > 0. From (|43t . 
we then get 

(4.5) WAT = I - N-^ + 0{N-^) mod 27rZ. 
From (14. 4p . we have 

g2«iv _ 2 (riV + 0{N-^)) e"" - 1 = 0, 

which implies 

(4.6) VN = logi2TN) + 0{N~'^). 

Since f{z,T) ~ 2arcsin(zr) and arcsin(— z) = — arcsin(z), (14. 5 1 and (|4.6p prove the 
lemma. □ 

We define the ^^-norm oi n — {ni, • • • , n^) e Z'^ by 

d 

(4.7) |n|,i=5]K|. 

i=i 

We also introduce the following notation. For n £ Z'^, we define B{n) £ 
B(Z2(Z'');C) and K{n) G B(C;;2(Z'^)) by 

(4.8) Bin)f = fin), 

(4.9) {K{n)c){m)^c6„,n. 
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Note that 

(4.10) P{n) = K{n)B{n) 
Let ro(fc, z) be defined by 

(4.11) ro{k,z) ^ {2n)-'^ I —— dx 



rpd h{x) — z 

Then the resolvent Rq{z) — {Hq — z)^^ is written as 

(4.12) (i?o(z)./)(m) = ^ ro(m-n,z)/(7i). 

Lemma 4.2. For any m,n and \z\ > d 

oo 

(4.13) B(m)Ro{z)K{n) = ^ Cs{m - n) z'"'^ , 

s—\7n—7i\ 

where Cs{k) is a constant satisfying 

(4.14) |cs(fc)|<c^^ s = 0,l,2,---. 
In particular, for \z\ > 2d 

(4.15) \\B{m)RQ{z)K{n)\\ < \z\-^-^"'-''^'\ 

Proof. Using {h{x) — z)^^ — — X)^o z^'^^^h{x)'' , we have for large \z\ 



oo 



(4.16) ro(fc,z) = ^z-^-ic,(fc), 



s=0 



Cs{k) ^ -(277)-" / h{xYe"'-''dx. 
Jt'' 

Since |/i(a;)| < d, we have 

(4.17) \cs{k)\<d'. 

Hence the series (|4.16p is absolutely convergent for \z\ > d. Note that h{x)'' is a 
sum of terms of the form 

(cosxi)"^ • • • {cosxd)"'^ , < ai + ■ ■ ■ + ad < s, < aj < s. 

If s < \k\ii , we have aj < \kj\ for some j, which implies 

/ h{xye'''-'=dx^O, if s<|fc|;i. 

Then we have 

(4.18) cs{k)^0, if s<|fc|,i. 
By (|4l^ and filB . we have 

OO 

B{m)Ro{z)K{n) = ^ z"*"ics(m - n). 

Then lemma then follows from (I4.17P and (|4.18p . □ 
Lemma 4.3. For any m, n, there exists a constant C'mn such that if\z\ > \\H\\ + I, 
\\B{m)R{z)k{n)\\ < C„„,(l + |zri-l™-"l'\ 
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Proof. Let p = \m — n\ii . By tlie perturbation expansion, we have 

- Roiz) - Ro{z)VRo{z) + ■■■ + {~lYRo{z)V ■ ■ ■ VRo{z) + 0{z-P-^). 

Multiply this equality by B{m) and K{n). Then by Lemma 4.2, the first term 
decays like 0{\z\~^~^) . Next we look at the term 

B{m)R(i{z)V ■ ■ ■ VRo{z)K{n), 

consisting of k numbers of V and fc + 1 numbers of Ro{z), where 1 < k < p. It is 
rewritten as a finite linear combination of terms 

(4.19) B{m)Ro{z)P{A^'>)Ra{z)P{A^^) ■ ■ ■ P{r^^^)Rn{z)K{n), 
since V — X]|r|<c ^(^)^('') for some c > 0. We put 

ei = |m-r(i)|,i, = jr^D - r'^) , • • • , = [r^ - . 
By (|4.10p . (|4.19l) is written as the product 

B{m)Ro{z)k{r'^^^) ■ B{A^'>)Ro{z)K{r^^^) B{A''^)Ro{z)K{r^'^^). 

By (I4.15p . this decays like |2;|-('''+i+^i+ - +<^fc+i). Since \m- = p, we have 

ei H h Efc+i > p. 

Taking notice of fc + 1 + ei + • • • + Cfc+i > 2 + p, we have proven the lemma. □ 

We can now solve the inverse problem for H. 

Theorem 4.4. Suppose V is compactly supported. Then from the scattering am- 
plitude A(X) for all X € (0,d) \ (Z U ap{H)), one can reconstruct V . 

Proof. Let A{\; 6, 9) be the integral kernel of A{X). Let \/A = fc, and put 

Bik; 0, 9') = {J{ke)J{k9'))-' A{k^;9, 9'). 

In view of (|3.18p . p.22p and (13.231) . we can rewrite it as 

(4.20) B{k- e, 9') = Bo{k; 9, 9') ~ Bi{k; 9, 9'), 

(4.21) Bo{k-9,9')= e-™-(="('='')-="('='''»i>(n), 

(4.22) Bi{k-9,9') = - e-"-"('=^)F(n) (^(fc2+iO)F(?(°)(fc,0')) 

(4.23) ^(o)(fc,6l') = fe'"''^^'''''^! 

Let 9j ^0, 9'j^O {] = !,■■■ , d), and put C(z, 9) = (/(z, 9i), • • • , /(z, 9d)), where 
f{z,T) is defined in Lemma 4.1. Then Bo{k;9,9') and Bi{k;9,9') have analytic 
continuations Bq{z;9,9') and Bi{z;9,9'), which are defined with k replaced by z 
in the upper-half plane and x{k9) by C{z,9). We put 

S{n) = ni + ■ ■ ■ + rid, n G Z'*. 
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We now take 9j > and 9'j < 0, {j = I,-- - ,d). Then by Lemma 4.1, as 

z ^ N + i ^ oo, 

(4.24) e-'"-*^'^'^) = e-»^(")^(47V2)^(") n {0j)^"' (1 + 0{N-^)), 

j=i 

(4.25) e'"<^'-'''^ = g^s(«)7r(4^2)S(«) ^ ^^^^2", _^ O(iV-i)). 

The recontsruction procedure for V{n) goes inductively with respect to S{n). 
Let us assume that 

(4.26) suppy C {(m,--- ,nrf); |n,|<M}. 

We use (|4.2ip . (|4.24p and (14.25^ to compute the asymptotic expansion of Bo{zn; 9, 9'). 
zn — N + i, a,s N oo. Then the largest contribution arises from the term for 
which S{n) is the largest, i.e. n = (M, • • • , M). Therefore, we have 

(4.27) BoizN ; 9, 9') ^ {2Nf'''' n (0,0;)'*'i>(n(*^)), 

where n^*^) (M, • • • ,M). Since 

(4.28) \\R{{N + if)\\^0{N^^), 
using ([422]) . (|424|) and (|4?25|) . we see that 

(4.29) Bi{zN ; 0, ^^') = 0{N^-^^''-^). 

By (|4.27p and (|4.29p . we can compute y(n*^*^^) from the asymptotic expansion of 
B{zn;9,9'). 

Assume that we have computed V{n) for S{n) > p. Then 
Bo{zn;9,9')- J2 e-*"-(^(^''')-';("'«'»y(n) 

S{n)>p 

(4.30) ^ 

The image of the map 

gd-i ^ gd-i 3 ^ ^ . . . ^ 

contains an open set in R'', where 

5-^-1 = {61 e S'^-^ ; 6lj > 0, Vj}, S^-^ = {0 ^ 5"^"^ ; 9j < 0, Vj}. 

Therefore, one can compute V{n) for S{n) = p from (|4.30p . 

We show that Bi{zn',9,9') = 0{N'^p^'^) up to terms which are already known. 
We rewrite Bi{zn;9,9') as 

B.Izn; 9, e') = e-^"<(^~'«)e"-^(""^^') • V{m)V{n) ■ B{m)R{zl)k{n). 

m.n 

We split this sum into 4 parts: 

E + E + E + E =--h + h+h+h. 

S(m),S{n)>p S{m)>p,S{n)<p S{m)<p,S{n)>p S{m),S{n)<p 



DISCRETE SCHRODINGER OPERATORS 



17 



Note that by (I4.24|) and (|4.25p . we have 

Then by and (|i3T]) . we have 

(4.32) /4 = 0(7V^P~2). 

By Lemma 4.3 and (|i3T]) . I3 = o(7v2(s(m)+s(«)-i-|m-n|,i)-)^ gj^^^g 

S{n) — S{m) ~ ^^(%' ^ ™j) < I^M — , 

we have 

2S'(m) + 2S{n) - 2\m ~ n\ii - 2 

= 45(to) + 2{S{n) - S'(m) - |m - ) - 2 < 4p - 2, 

which proves 

(4.33) I3 = 0(7V'*f-2). 
Similarly, we can prove 

(4.34) I2 = 0{N'^P-^). 
We finally observe Ii . We put 

S(ti)<p S(n)>p 
i?>p = + V>p, Rypiz) = {Hyp - z)-\ 

By the resolvent equation, Ii is split into 2 parts 

/i= E e-™-^(""'*)e""'^(""'^')-F(m)F(n)-B(TO)i?>p(z^)i?(n) 

S(m),S(n)>p 
S(m),S{n)>p 

The 1st term of the right-hand side is a known term, since we have already recon- 
structed V{n) for S(n) > p. The 2nd term is a linear combination of terms 

(4.35) ^-^m■c{z«.e)^^n■ciz^,e') . B {m)R>p{z%)P {k)R{zl)K (n) , 

where S{k) < p. By Lemma 4.3, it decays like 0(iv2(^(™)+s(n)-|™-fc|,i-|"-fc|,i-2))^ 
Using 

S{m) - |m - k\ii = E("^J - ^j) - E l™J " '^jI + E 
j j J 

<5(fc)<p, 

we see that (|I35)) decays like 0{N'^p-'^). Therefore, we have 

(4.36) h = OiN^P-^), 

up to known terms. By virtute of (I4.32p ^ (|4.36p . we have completed the proof of 
the theorem. □ 
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5. Estimates of the free resolvent near the critical values 

The purpose of this section is to derive estimates of the resolvent Ra{z) = {Hq — 
z)~^ in weighted Hilbert spaces. Equivalently, we consider the operator norm on 
i^(Z'') of qRo{z)q, where q is the operator of multiplication by g e £°°{Z'^) : 
(qf){n) = q{n)f{n). In particular, ps is the operator defined by 

p,H = (l + |n|2)-V2^ seR. 



We put 



Di = {z e C; < |z| < 1}, 
Ad = p(i?o) = C\[0,d], 



(5.1) A(z) = i(^2-z-i 

For w = re*^ with < d < 2n, we take the branch ^/w = y^e*^/^. For a, 6 € C, we 
put 

I{a, b) ^ {ta+{l-t)b;0<t< 1}. 

Lemma 5.1. X{z) is a conformal map from Di onto Ai, and its inverse is given 
by 



(5.2) z(A) = 2A-1-2VA(A-1). 

Proof. Since the map 

{0 < \z\ < 1} 3 z ^ l ( z + -] e C \ [-1, 1] 



2 

is conformal, so is A(2;) from Di onto Ai. By solving the equation + (4A — 
2)z + 1 = 0, we have the inverse map z = 2A — 1 ± 2-y/A(A — 1). For A > 1, 
|2A - 1 - 2^A(A - 1)1 < 1. Therefore, we obtain jSj]). □ 

The following lemma is proved by Lemma 2.4 and Theorem 2.6 (2). Note that 
Hq has no eigenvalues. 

Lemma 5.2. For s > 1/2, the operator-valued function ps Ro{X) Ps is analytic with 
respect to X G Ad, and has continuous boundary values when X approaches E ± iO, 
Ee (0,d)\Z. 

We study estimates for Ro{E±iO) when E e (t{Ho) = [0, d] is close to {0, 1, • • • , d}, 
the set of critical values of h{x). Let us begin with the case d — 1. 

Lemma 5.3. Assume d — 1. 

(1) roln, z) defined by ^-11^ has the following representation 

^"^"''^^^^ = ^^^' VA(.)(A(.)-1) ' K-)eZxD, 

Moreover, ro{n, X{z)) has a meromorphic continuation from Di into C. 

(2) Let II • Whs be the Hilbert- Schmidt norm on ^^(Z), and take qj — {qj{n))ne'L € 
i!^(Z), i = 1, 2. Then we have for A G Ai 

\\nR(\\n\\ ^ llglll^^(Z)l|g2||£^(Z) 

UR.{X)q,\\us<^^^j—=^. 
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M(A,Ai)=(l + 2 max ) ( 1 

A^(A,Ai) 



(3) Let u{v) = v{v — 1), and take A, A G Ai such that /(A, Ai) C Ai. We put 

1 ^ 1 \ 
'veilTxi)^^^"'^) y ' \u{\)\ |u(Ai)|y' 
1 1 
u(A) u(Ai) 

Then for any < a < 1, we have the following pointwise Holder estimate 

(5.3) |ro(n,A)~ro(n,Ai)| < |A - Aip (1 + |n|)"Af (A, Ai)"7V(A, AO^"", 
and the following Holder estimate of the Hilbert- Schmidt norm 

(5.4) WMi (RoW - Ro{^i))Paq2\\HS < |A - AirC„(A, Ai)||gi||,2(z)||g2||£2(z), 

(5.5) a(A,Ai) - M(A,Ai)"iV(A,Ai)i-". 
In particular, there exists a constant C'^ > such that 

WPaqi (RoW - -Ro(Al)) Paq2\\HS 

(5-6) C"|A-Ai|" 

- |A(A-l)Ai(Ai-l)|(i+3")/2 \\l^MW2hHz), 

2/|A|,|Ai| <2, anrf/(A,Ai)c Ai. 

Proof. To prove (1), we first note by residue calculus 

rQ[n, Xiz)) ~ — I — ^ — — / — -— ■ — — . 

TT Jo i — cosx— 2A(z) TTJ (w — z)(w — i/z) z—l/z 



By ([52]), we have i/z = 2A - i + 2^A(A - i). Hence z-l/z = -4^A(A - f), 
which proves (f). 

Using (|4.f 2|) and (f ), we obtain for z G Di 

||gii?0(A(z)) 92111,5= E |;,(J)(;,(^)„i)| lgl(")ng2(HP 



|A(z)(A(z)-l)|' 

where ||qj||2 = Y.ne7. \qj ■ This proves (2). 

We prove (3). Let Ai, A G C+ and z = z(A), zi = z(Ai). Then (1) and (|5?2|) give 

zH 4"! z'-l-zH , y,\ u{X)-u[\^) 

ro n,A) -ro(n,Ai) = -— + — — = — + z[ — . 

u[\) u[\i) u[\) u{X)u{Xi) 

Since z, zi G Di, we have 

|n|-l \n\-l 

-zl"l| = l(z-zi) J2 ^'"'"'"'^il < E 1| - l^ll^-^il- 

j=0 j=0 

Moreover, we have by using ([57 



KA)-«(Ai)| < ( max \u'{iy)\ ] |A-Ai|, 

, l^Gi (A.Ai ) 



|^(A)-z(Ai)| <2 i+ max K(i.)| |A-Ai|. 
i/e/(A,Ai) 
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The above inequalities yield 

|ro(n,A)-ro(n,Ai)| < |A-Ai| ( + max U {,,)\ (2\n\ + --\-r)) . 

\\u{X)\ iye/(A,Ai) V |u(Ai)|// 

Interchanging A and Ai, and adding the resulting inequalities, we have 

(5.7) |ro(n,A)-ro(n,Ai)| < |A - Ai |(1 + |n|)M(A, Ai). 
We also have by (1) 

(5.8) |ro(n,A)-ro(n,Ai)| <iV(A,Ai). 

By dSJ]) and we obtain ((O)) . 

Using and (fO)) . we have for z e Di 

WPaQl (^o(A) - i?o(Al)) Paq2\\HS 

= X! \'^o{n- m,\) - ro{n- m,Xi)\'^\pa{n)qi{n)\'^\pa{m)q2{rn)\'^ 

< |A-Aip"(l + |n-m|)2"C„(A,Ai)2|p„(n)gi(n)np„(m)g2(m)p 

< |A-Aira(A,Ai)2 J2 l9iHn$2(m)p, 

which proves (|5.4p . The inequality (|5.6|l follows easily from (|5.5p . □ 

We study the case d = 2. 

Lemma 5.4. Let d — 2 and q{n) — qi{ni)q2{n2), where qj £ £^(Z) and n = 
(711,712) S Z^. T/ien t/iere exists a constant C > smc/i that 

(5.9) ||?i?o(A)Ql| <C||gi||2.(2)||g2L'.(z)|fog(A(A-l)(A-2))|, 

/or aZ/Ae A2n{|A| < 3}. 

Proof. We prove the lemma by passing it on the torus. The idea consists 
in reducing it to the 1-dimensional case, regarding the remaining variable as a 
parameter. We put 

g,(x,) = (2vr)-i/2 ^ $,(",) e»^-^ 

and define the convolution operator qj* by 

{q, * f) (x) = / q,{x, - y,)f{y)dy„ f E L\T'), 
Jo 

where y = (2/1,2:2) if j = 1, y = (2:1,2/2) if j = 2. We put 

fi = fj.{X,X2) = A - h{x2), h{x2) = ^ (1 ~ COSX2) , 

and define the f-dimensional operator Ai{p) with parameter /i by 

yli(/i) = qi* {h{xi) - qi* ^qi* {Hq - A)"^ qi * . 

Take /,/' e L2(T2), and let /, /' e £^(22) be their Fourier coefficients. We are 
going to estimate 

C2(A) := {qi 92 (^^0 - A)"^qi % /, ]')t:\'L^) 

= {q\ *q2* {Ho - Xy^qi *q2* /,/')i2/T2) ■ 
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Letting 

5 = 92 * /, 5' = 52 * /', 

we have 

C2W = J ((yli(^)g)(-,X2),.9'(-,X2))^2(Ti) '^^2- 

By Lemma 5.3 (2), we obtain 

ll5(-,a.2)||L^(Ti)||<7'(-,a;2)|U=(Ti) 



Since 

II (92 * /) (•,a:2)||L2(Ti) < ||g2|U2(Ti)ll/IU2(T2), 

which follows from a simple application of Cauchy-Schwarz inequality, we have 

|C2(A)| < |kl||i.(Ti)ll92||i2(Ti)ll/l|L^(T^)ll/'l|L^(T^)^2(A), 



D2{X) 



dX2 



/o |/i(A,a:2)(M(A,X2)-l)|l/2- 

Then the lemma is proved if we show for |A| < 4 

(5.10) i:>2(A) < 8 {K{\) + K(X - 1) + K{\ - 2)) 

3^/2^ 



X(A) = V2|^2 + log:|^l 



To prove (|5.10p . we let 



1 



J(A) = / \{s-X)s{l-s)\-^/^ds, 



and first derive 

(5.11) i:'2(A) < 2J(A) +2J(A- 1). 

In fact, by the change of variable s = (1 — cosa;2)/2, we have 

D2{X) =2 f \{s- X){s + 1 - A)s(l - s)\-^/^ds. 
Jo 

Using the inequality 

< 



(5.12) ^ 



a a — 1 



|a(a- l)|i/2 

with a = s + 1 — A, we obtain (|5.1ip . 
In order to compute J(A), we put 



1/2 ^ 



i|i/2 |a_l|i/2 



"1/2 

Jo(A)= / \s{s - X)\-'/^ds. 



Then we have 

^1/2 , d 

J{X)<2 \ +2 / ^ = 2Jo(A)+2Jo(A-l). 

Jo |s(s-A)|2 7i/2 |(l-s)(s-A)|i/^ 

This, combined with (j5.11l) . implies 

(5.13) D2iX) < 4 (Jo(A) + 2Jo(A - 1) + MX - 2)) . 
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Let US first consider the case |A| > 1. Estimating as in (|5.12|) . we have for 
< s < 1/2 



|s(s- A)|l/2 - |A|l/2 \^|s|l/2 |s_ A|l/2y - sl/2 (l_s)l/2- 

Hence for |A| > 1 



1/2 



1 



(5.14) Jo(A)<y^ ^^+(13^;'^^ -2. 

Next we consider the case |A| < 1 and let A = /i + iu. If > 4|/i|, we put 
s = \i^\t, E = l/(2|z^|), /lo = ^i|\v\. Note that y/2\v\ > |A| > and if t > 1/2 then 
t- ^J.a> t/2, hence |(t - ^o)^ + ll^^^ > (t/2)i/2. We now compute 

.o(A) " 

(5.15) 



Vi|(<-Aio)2 + l|i/4 

dt dt 



< 



Vi|(<-Mo)2 + l|V4 71/2 Vt|(i-Mo)2 + l|i/4 

f'^'^ + V2 r ^ = V2 + ^/21og(2i^;)<2^/2 + ^/21og^. 
Jo Ji/2 r |A| 



'1/2 

If < 4|^|, we let s = \^\t, R = 1/(2|^|), a{p) = ^/|^|, j/q = v/\^i\. We then 
obtain 

Jo (A) 



x/t|(f-(7(/.))2+z.o|i/4 



(5.16) 



< 



Vt\{t-a{^,)Y+vl\^l^ J^/2^^t\{t-a{^i)Y + vl\y^ 
1/2 dt /-^ d< 



Vt|t-l|l/2 7l/2i-l 

3v/2 



< 2 + log2i? < 2 + log- . 

|A| 

In view of (|5.14l) . (|5.15p and (15.161) . we have 

Jo(A) < 72(^2 + log 

which, together with (|5.13l) . proves (|5.10p . □ 

Finally we consider the case d > 3. 
Lemma 5.5. (1) Let d > 3 and 

Qs{ni, ■■■ ,nd) ^ qi{ni)q2{n2)ps{n'), 

where n' = {123, • • • ,nrf) G Z'^"^ 91,92 G i'^i'Z) and Ps{n') = (1 + |n'|2)-''/2 with 
< s < 1. Then there exists a constant Cs > such that the following estimate 
holds: 

(5.17) \\Qs%{X)Qs\\ < Cs\m%iz)\M%iz) 

for an A e Ad n {A e C ; |A| < 2d}. 
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(2) Moreover, let A, Ai G C± n {A e C ; |A| < 2d} and let g > be small enough 
and 

Qs.gini, ■■■ ,nd) = Pgini)Qs+2g{ni, ■ ■ ■ ,nd), n G Z"*. 
Then there exists a constant C's.g > such that the following estimate holds true: 

(5.18) WQ^.g (i?o(A) - i?,o(Ai)) QsJ < |A - Xi\'CsJqi\\%^^^\\q2\\%iz), 
Proof. (1) The proof is similar to the one for the previous lemma. Let 

^2(/^) = Q'l * 92 * (-f^o — A)~^gi * 52* = gi * 92 * {hi + /12 — * 12*, 

1 

hj = h{xj) — ~ l^ ~ '"('^' x') = X — ^2 ^j' 

i=3 

where x' = {xs, ■■■ , xa). For /, /' e L^{T'^), we put gs= Ps* /, g's = Ps* /'• Then 
we have 

C(A) := (gi ?2 (^0 - A)-igi % pf, pf) 

= {qi * 92 * {Ho - X)~^qi * q2* Ps* f,Ps* /') 

= / {A2{lJ')gs{-,x'),g's{-,x'))L2(^T^)dx'. 

J'rd-2 

Lemma 5.4 then implies 

\c{x)\ < c\M%^^02\\%^z) 

X / |log(/u(/U-l)(/U-2)|||5,(-,a;')||L2(T2)||5s(-'a;')l|L2(T2)rfa;', 
where C is a constant independent of A e A^. We now put 

Ds{X) = [ |log(/i(A,x')(M(A,x') - l)(/i(A,x') - 2)|i|(p, * /)(•, x')|ii.(T2)dx'. 

Lemma 5.5 will then be proved if we show the existence of a constant Cg independent 
of A e Ad n {|A| <2d} such that 

(5.19) D,{X)<C,\\f\\l,^^,y 
We define the set SP by 

SP = {(a;3,--- ,Xd); Xj = or tt, j = 3,--- ,d}. 

This is the set of singular points for ^{X,x'), since 'Vx'IJ'{X,x') = if and only if 
x' G SP. We label the points in SP by p<-^\- ■ ■ N = 2'^-'^: 

SP = {pW,...,pW}. 

For a sufficiently small e > 0, we put 

= {x' e T'^-^ ; \x' - p^-?') I < e}, 1 < j < A^, 

T(0)(A)=T'*-2\(uf=^T«), 

and let 

Ei^\X) = [ |log(Ai(A,x')(M(A,a;') - l)(/x(A,x') - 2)|||(p, * f){;x')\\l.^T..^dx' , 
E^J'^\X) = [ \log{^,{X,x') - fc)|||(p, * f){;x')\\l.^^.^dx'. 
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Then we have 

N N 2 

(5.20) D,{\) = Y.E(^\\) < ^5]i?p''=)(A). 

j=0 i=0 fc=0 

We shall make use of the following version of Heinz' inequality. For the proof, 
see p. 232 of [6 . 

Proposition 5.6. Let A, B be self-adjoint operators satisfying 

A>1, B>0, D{A)(zD{B), \\BA-^\\<1. 
Then for any Q < 9 < \, we have 

D{A^) C DiB"), \\B'^A-^\\ <1. 

Let us estimate Ei'''''\x). We take 6 = -, and define self-adjoint operators 

A and B in L^{T''-^) by 

Af(n') = Co(l + |nf)("-2)/27(n'), 

(Bf) {x') = x^^\x')\ \og{fi{X,x') - k)\'^'f{x'), 

where x^''H^') is the characteristic function of the set T^^^ and Co is a constant to 
be determined later. We compute the Hilbert-Schmidt norm of BA~^. Let k{x') be 
the inverse Fourier image of + Then, up to a multiplicative constant, 

BA^^ has the integral kernel 

(5.21) Kix',y') = Cox'^^\x')\ log(/x(A,x') - k)\'^'k{x' - y'). 
One can show that for any a > 1, 

(5.22) sup / |log(^(A,x') - fc)|"da;' < cx). 

AGAdn{|A|<2£;} JTj 

In fact, if i — 0, Va;'yu(A, x') 7^ on Tq. Then we can take X]j=3 cosx^ as a new 
variable to compute (|5.22l) . The case j ^ is dealt with as follows. Suppose, for 
example, p^^-* = (0, • • • ,0). By the Morse Lemma, we can introduce new variables 
3 < j < d, around p'^^-' so that J2j=3 cosxj — (d — 2) — J2j=3 v]- One can then 
prove (|5.22l) by an elementary computation. The other cases are treated similarly. 
On the other hand, by Parseval's formula 

/ |fc(x')|'dx' = y (1 + |nf )-("-2) ^ 

Therefore BA^^ is of Hilbert-Schmidt type, in particular, bounded. By choosing Co 
small enough, we have ||i3v4^^|| < 1. Then by Proposition 5.6, B^A^^ is bounded 
on L^(T'*~^), which implies that 

||£;PHA)/|| <C.||/||i.(T<<), 

This proves ([5l^ . 

The proof of (2) repeats the arguments from the proof of (1). □ 
As a consequence of the above lemma, we show the following theorem. 
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Theorem 5.7. Let d > 3, and = Hq + 'yV, where V is a complex-valued 
potential such that V{n) — 0(|n|^'*), s > 2, as \n\ — > oo, 7 being a complex 
parameter. Then there exists a constant 6 > such that has no eigenvalues 
when I7I < 5. 

Proof. We put Q{ni,--- ,nd) ^ qi{ni)q2{n2)p{n'), where 
gi(ni) = + 52(^2) = (l + |n2^-('+^)/^ P(n') = (l + |nf)-^/^ 

If there exists S e R and / e ^^(Z'') such that {H^ + 7^ - E)f^ 0, we have 

Choosing e > smaU enough, we have Q^^ V Q^^ E B(i?^(Z'')). Using Lemma 5.5 
and taking 

we obtain Theorem 5.7. □ 

6. Traces formulas 

6.1. Predholm determinant. We shall discuss in an abstract framework in this 
section. Let "Hhe a. Hilbert space endowed with inner product ( , ) and norm || • ||. 
Let Ci be the set of all trace class operators on H equipped with the trace norm 
II • llcj. Recall that for K € Ci and z G C, the following formula holds: 

(6.1) dot (I - zK) =cxp(^- Tt{K{1- sKy'^)ds^ 

(see e.g. [5]). As is well-known if A e B('H;'H) and B G Ci, we have 

(6.2) Tr(AB) = Tr(BA), 

(6.3) det(/ + AB) = dct(/ + BA). 

Suppose we are given an operator H = Hq + V on H satisfying the following 
conditions: 

(B-1) Hq is bounded self-adjoint. 
(B-2) V is self-adjoint and trace class. 

We put 

Ro{X) ^ (Ho - X)-\ XeA:=piHo), 

and define D{X) by 

(6.4) D(A) ==det(/ + yi?o(A)), A e A. 
Lemma 6.1. (1) D{X) is analytic in A. Moreover 

(6.5) L'(A) = 1 + 0(1/A) as |A| ^ 00, 

(6.6) iogD{X)^-Y,^-^TT {VRo{X)r , 

n=l 

where the right-hand side is absolutely convergent for \X\ > tq, tq > being a 
sufficiently large constant. 

(2) The set {A £ A; D{X) — 0} is finite and coincides with ad{H). 
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Proof. Letting Ev{t) be the spectral decomposition of V , we define 

/•OO /"OO 



There exists ro > such that ||i?o(A)|| < C/|A| for |A| > ro, which imphes 
(6.7) |Tr (T/i?o(A)) I < C/|A| for |A| > tq. 

By (|6.2p . taking ro large enough, we have for |A| > tq 
Tr (Fi?o(A)(l + st/i?o(A))-i) 
= Tr (yi/2i?o(A)(l + 5yi?o(A))-i|y|i/2' 



= ^(-s)"Tr (yi/2i?o(A) VRoiX) ■ ■ ■ VRoiX) \V\'/^) 

OO 

= ^(-.rTr(yi?o(A))"+\ 

We then have (|6.5p and (|6.6p by (|6.ip . For A ^ cr{Ho), the eigenvalue problem 
(iJ — A)m = is equivalent to (/ + {Hq — X)~-^V)u = 0, which has a non-trivial 
solution if and only if det(/ + {Hq - Xy^V) = 0. This proves (2) by dO]). □ 

Lemma 6.2. The following identity holds: 

(6.8) logDiX)^-J2^(\) , F„=Tr(i/«-i7o"), n>l, 



where the right-hand side is unifomly convergent on {|A| > rg} for Tq > large 
enough. In particular, 

(6.9) i^i = Tr(y), F2^Tr{2VHo + V^). 

Proof. Let R{X) = {H — A)^^. Take rg > large enough. Then for |A| > ro, we 
have by the resolvent equation 



(6.10) i?(A) - i?o(A) = ^(-l)"i?o(A) VRoiX) ■ ■ ■ VRo{X) . 



n=l 



Let F(A) = logi:»(A). Since jtRo{X) = i?o(A)^ we have by (HH) and 



n— 1 

= Tr(i?(A)-i?o(A)). 

Using the equation 

OO 

i?(A) = - ^ A-"-ii/", 

ri=0 

we obtain 

OO 

F'{X) = ^A-"-iTr(i?"-H^). 

11=0 

In view of ([6?5|) . we get ([6?8l) . □ 




Figure 1. The spectral shift function 5(A) and five eigenvalues h, - ■ ■ ,1^ for H 



6.2. Spectral shift function and trace formula. Let H = Hq + V satisfy (B-1), 
(B-2). Then there exists a function ^(A) such that the foUowing equahty 

(6.11) Tr ifiH) - /(i/o)) - / e(A)d/(A) 

JR 

holds, where / is an arbitrary function from some suitable class. If / is absolutely 
continuous, then df can be replaced by f'{X)dX. We call (|6.1ip a trace formula, 
and ^(A) the spectral shift function for the pair H,Hq. A typical form of ^(A) for 
our case of discrete Schrodinger operator H in §2 is drawn in Figure 1. 
Let us recall the basic properties of ^(A) (see [7], [31] )■ 

(1) The following identity holds: 

\ogD{X) = J^&-dt, AeC+, 

where -D(A) is the perturbation determinant defined by ^6.4^ , and the branch of 
\ogD{X) is chosen so that log-D(A) — o(l) as |A| oo, and ^{t) G _L^(R). We 
have 

^(A) = lim — argZ?(A + ie), a.e. A e R, 

£-i-+0 TT 

where the limit in the right-hand side exists for a.e. A G R The support o/^(A) is 
equal to (t{H) and 

I m)\d\<\\v\\c,, 

Jr 
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(2) Its relation to the S-matrix is 

(6.12) det5(A) e-^'r^CC^), for a.e. A e aac(Ho). 

(3) If H have > negative eigenvalues and N-^- > positive eigenvalues, then 

- 7V_ < ^(A) < iV+, fora.e.XeR. 

(4) Suppose Hq has no eigenvalues in the interval (a, &) C R. Assume that Xq G 
(a, 6) is an isolated eigenvalue of finite multiplicity dg of H . Then ^(A) takes an 
integer value n_ (n^^) on the interval (a. Aq) (on the interval (Xo,b)). Moreover, 
we have 

(6.13) C(Ao + 0) - e(Ao - 0) = -do. 

(5) IfV>0, then ^(A) > for all A e R. 

(6) IfV<0, then ^(A) < for all A e R. 

(7) If the perturbation V has rank N < oo, then ~N < ^(A) < N for all A e R. 

As will be shown in the following lemma, Fn/n, the Taylor coefficients of — log\D(A) 
around A = oo are equal to the moments of the spectral shift function ^(A). The 
first two terms were computed in Lemma 6.2. To compute the terms for n > 3, we 
impose the following assumption. 

(B-3) There exist unitary operators Sj (1 < j < d) such that 
1 

(6.14) Ha = --{s + S*), SjS^ = S^Sj, Vz,j, 

(6.15) Tr{S^VP)=0, Vj = l,..,d, fc ^ 0, p>l. 

Note that by KT^i and (|6T5l) . we have 

(6.16) Tr{S''VP)=0, Tr (^^(5'*)'=) = 0, fc > 1, p>l, 

(6.17) Tr(S'°(S'*)Vf) =0, a b, p>l. 
Lemma 6.3. Let H = Ha + V satisfy (B-1), (B-2) and (B-3). Then 

(6.18) F„ = Tr (i/" - H(i) = n [ e(A)A"-MA, n>l. 
In particular, letting t = —1/4, and 

d 

(6.19) I\V ^tY^(SjVS* + S*VSj), 

i=l 

we have 

(6.20) Fi=Tr(V^), F2^TT(y^), F3 ^ Tr {V^ + edr'^V) , 

(6.21) F4 Tr (V^ + SdT^V^ + 2t{AV)V) , 

(6.22) F5 = Tr (V^ + 30d(2d - 1)tV + lOdrV'"^ + 5t{AV)V'^) . 
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Proof. By taking /(A) = A'' in (16. lip , we get (|6.18p . In Lemma 6.2, we have 
proven that Fi — Tr (V"). F2 and F3 are computed by the use of (|6.15|) as follows: 

F2 = Tr {{Hq + Vf - Hi) = Tr {2H(,V + V^) = Tr [V^), 

F3 ^ Tr {{Ho + Vf - Hi) = Tr {SH'lV + 3HoV'^ + V^) = Tr {SH^V + V^), 

Tr {HlV) = r^Tr ((5^ + 255* + S*^)V) = 2T^Tr {SS*V) 

^ 2T^Tr {{d + SiS*)V) = 2dr2Tr {V). 

■Mi 

To calculate F4, we first compute 

F4 = Tr {{Ho + V)' - H^)) 

= Tr {V^ + AH^V + AHoV^ + AHIV^ + 2HoVHoV). 

Due to (|6.16p . we have 

Tr {HlV) = 0, Tr {HqV^) = 0, Tr {H^V'^) = 2dT^TT {V^). 

Using Tr {SVSV) = 0, we get 

Tr {HoVHoV) = r^Tr ((5 + S*)V{S + S*)V) 

= T^Tr {{SVS* + S*VS)V) = rTr {{AV)V,), 

since Tr{Y.^^^S,VS*V) =0. 
Finally we compute F^. Firstly, 

F5 = Tr {{Ho + V)^ - Hi) 

= Tr {V^ + 5HlV + 5HoV^ + 5HIV^ 

+ 5HlV^ + bHlVHoV + bHoVHoV^). 

By (j6T6| . we have Tr(i?oV4) = 0, Tr(i?5^y2) ^ q, and Tr (i/^Vi/o^^) = 0. We 
then have 

Tr (HoVHoV^) ^ r^Tr ((5 + S*)V{S + S*)V^) 
= r^Tr ((5^5* + 5*^5)^2) 

= rTr {{AV)V^) ^ Tr (hIv) = r^Tr ((5 + S*)^V) 
= 6t^Ti- (5^5* V) = 6T4d(2d - l)Tr (V^). □ 

The above lemma enables us to estimate the eigenvalues in terms of the spectral 
shift function. 

Theorem 6.4. Let H = Ho + V satisfy (B-1), (B-2) and (B-3). Assume that 
(j{Ho) = [a, /3], and put 

(6.23) Er,^n [ C(A)A"-idA. 

iR\[a,/3] 

Let nij he the multiplicity of Xj G (Jd{H). Then we have for any n > 

(6.24) = 
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(1) IfV>0, then ad{H) C oo) and 

(6.25) A,"^,<Tr(y) 

(6.26) A|™,<Tr(F3 + |^y). 

(2) IfV <0, then ad{H) C (-oo,a) and 

(6.27) ^ A,m, >Tr(F), 

(6.28) 5: A|™,>Tr(F-V|^F). 

Proof. For small i > 0, we define the set 

Ot = (-00, a-t)\j{l3 + t, oo). 

If \.j G ad{H), we take e > small enough so that the interval Ij = {Xj — e, Xj + e) 
satisfies Ij n ad{H) = {Xj}. Then we have, by the property (j6.13p . 

{X)dX — -'■mj5{X ~ Xj)dX, on Ij 

(see Fig. 1). More precisely, see (3.28), (3.29) of [7 . We then obtain 

En{t):^n[ £_{X)X"-^dX = - [ e'(A)A"dA 

= E AJ'm, ^£;„(0) = Y A7toj =£;„, 

which proves (I6.24p . 

If F > 0, then Xj > and ^(A) > 0, which implies E2n-i < i^2n-i- Then 
and (jOni) give (|05| and ([g?^ . The proof for the case V < is similar. □ 

Remark 6.5 For our discrete Schrodinger operator discussed in sections 1, 2, 3, 
V = J2neZ'^ Vn{n)P{n) is trace class if 

E \vi^)\<^- 

The assumptions (B-1), (B-2), (B-3) arc then satisfied if we shift our Hamiltonian 
Ho in §2 by d/2. 

Remark 6.6 For the continuous model, it is well-known that H has no embedded 
eigenvalues for the short-range perturbation. The (non) existence of embedded 
eigenvalues for the discrete model is an interesting open question. 
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